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We present calculations of the zero-temperature Casimir interaction between two free- 
standing graphene sheets as well as between a graphene sheet and a substrate. Results 
are given for undoped graphene and for a set of doping levels covering the range of exper- 
imentally accessible values. We describe different approaches that can be used to derive 
the interaction. We point out both the predicted power law for the interaction and the 
actual distance dependence. 
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1. Introduction 

Graphene^I^Sl has a very special band structure, leading to unexpected re- 
sults. Gomez-Santos^ showed that retardation effects are negligible for un-doped 
graphene. This is one of the peculiarities caused by the band structure. This means 
that we can limit the calculations to van der Waals interactions which are easier 
to derive. Throughout the text we make the basic derivations more general so that 
they are valid for any type of two-dimensional (2D) sheets, like, e.g., 2D metallic 
sheets. The result becomes specific for graphene when we insert the graphene po- 
larizability, a(q,u). In terms of the polarizability the dielectric function is given 
by e (q, lo) = 1 + a (q, u) = 1 — v 2D (q) \ (q,ui), where v 2D (q) = 2ire 2 /q is the 2D 
Fourier transform of the coulomb potential and x (<?> w ) the density-density corre- 
lation function or polarization bubble. For undoped graphene, in a general point in 
the complex frequency plane, away from the real axis the density-density correlation 
function iJ^ 

2 

X(q,z) = -J7TZ 7= =, (i) 
16ft ^Jv 2 q 2 - z 2 

where v is the carrier velocity which is a constant in graphene (E = ±hvk), and 
g represents the degeneracy parameter with the value of four (a factor of two for 
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spin and a factor of two for the cone degeneracy.) For doped graphene the func- 
tion becomes much more complicated^^ We refer the reader to Ref.[9] where the 
expression on the imaginary axis can be found in terms of real valued functions of 
real valued variab les. Ot her calculations where graphene were modeled differently 
has been publihedEM 

Section 2 is devoted to two freestanding graphene sheets, undoped or doped. We 
present derivations of the interaction and numerical results. In section 3 we treat 
a graphene sheet above a substrate. We derive the interaction and present results 
for a gold substrate. In section 4 we discuss the predicted power law dependence of 
the interactions and compare to the actual outcome from our calculations. Finally, 
section 5 we devote to summary and conclusions. Part of his work has previously 
been reported in a letter®. 

2. Two Freestanding Graphene Sheets 

It is possibleD^m^but very difficult to take spatial dispersion into account when cal- 
culating Casimir interactions in systems with three-dimensional (3D) layers. When 
dealing with parallel 2D sheets, like in the case of two freestanding graphene sheets, 
spatial dispersion enters in a natural way and does not cause any complications. 
Thus these systems can be used to gain experience concerning possible effects of 
spatial dispersion in general. 

2.1. Derivations 

The derivation of the non-retarded Casimir interaction can be done in several dif- 
ferent ways. Here, we will present two different derivations; one in the language of 
many-body theory; one in terms of electromagnetic normal modes. We begin with 
the many-body treatment. 




V 22 



Fig. 1. Feynman diagrams for the correlation energy in the two graphene sheet system. The 
ellipses represent polarization bubbles and the dashed lines the interactions indicated in the figure. 
The numbers 1 and 2 refer to which sheet the electron belongs to. See Ref. 15 for details. 
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2.1.1. Many-body theory 



For a system of parallel 2D sheets the derivation in the language of many-body 
theory becomes very simple, especially when retardation effects can be neglected. 
In the present two -she et system the interaction energy is nothing but the inter-sheet 
correlation energy^ISl. Using diagrammatic perturbation theory the Feynman dia- 
grams representing the correlation energy are given in Fig.[TJ To get the inter-sheet 
contribution we can either subtract the intra-sheet part or subtract the result when 
the separation between the sheets goes to infinity (at that limit only the intra-sheet 
contribution remains.) To avoid unnecessary complications we assume that the two 
sheets are identical. If we have two doped graphene layers the doping concentration 
is the same in the two sheets. It is straight forward to extend the treatment to differ- 
ent sheets. Each ellipse represents a polarization bubble, x(q, w), and the number 
indicates in which sheet the process occurs. Each of the four diagrams represent an 
infinite scries of diagrams. Since the sheets are identical the diagrams in the second 
row give contributions identical to the diagrams in the first row. We can use the first 
row diagrams and multiply the result with a factor of two. The interaction energy 
per unit area can be written as 

oo 1 

—^2 J £ J d\j2[Diagl(q,iu;;\)+Diag2(q,iu>;\)}, (2) 



o 



where A is the coupling constant and the factor of two has been inserted. 

Now vu = v 2D (q, ui) and v\2 — exp (— qd) v 2D (q, w), respectively, where d is the 
distance between the sheets. The exponential factor in the second potential is the 
result from taking the 2D Fourier transform of the coulomb potential in a plane the 
distance d from the center of the potential. The interaction lines with double bars 
represent a series of terms, with zero, one, two . . . number of polarization bubbles. 
This can be expressed as 

v n (q, w) = v 2D (q) + v 2D (q) x (q, cj) v n (q, w) + exp (-qd) v 2D (q) x (q, w) t?ia (q, w) 
Ci2 (q, ui) = exp (-qd) v 2D (q) + v 2D (q) X (q, w) v 12 (q, u) 

+ exp (-qd) v 2D (q) X (q, w) Uii (q, u) , (3) 

where we have closed the two infinite series. This system of equations can be solved 
and the result is 

- , , v 2D (g) [l + a(q,^)(l-exp(-2 g d))] 

Vu (q, lo ) = j 

[1 + a (q, w)] -exp(-2qd)a 2 (q, ui) 

~ , s v 2D (g) exp(-gd) 

Vi2(q,U)) = ^ • ( 4 ) 

[l + a(q,a;)] - exp (—2qd) a 2 (q, u) 

The square brackets in Eq. ([2} are [SuVnx 2 + C12W12X 2 — (v 2D /e) v 2D x 2 ], where 
the last term comes from the subtraction of the intra-band correlation energy. Each 
factor of e 2 appearing implicitly in the expression should be multiplied by the 
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coupling constant. Performing the integration over coupling constant gives 

(5) 



- e- 2 " d 



a (q, iui) 



1 + a (q, iuj) 



This is the result from diagrammatic perturbation theory within the random phase 
approximation (RPA). In next section we derive the same thing using normal modes. 

2.1.2. Normal-mode derivation 

An electromagnetic normal mode^^ is a solution to Maxwell's equations in absence 
of external perturbations. At zero temperature the interaction energy of the system 
is the change in the total zero-point energy of all the normal modes when interaction 
is turned on. In a planar system like the ones we treat in this work the modes are 
characterized by the 2D wave vector q in the plane of the sheets and substrates. 
Within this formalism the interaction energy is obtained as 

oo 

E c (d) = hj£^ j gln{/ q M}, (6) 



(i 



where / q (w q ) = is the condition for electromagnetic normal modes. One arrives 
at this expression after using an extension of the so-called argument principle and 
deforming the integration path in the complex frequency plane. 

We will now go through one possible way to find the normal modes and their 
zero-point energies. Let us assume that we have an induced carrier distribution, 
pi (q, ui), in sheet 1. This gives rise to the potential v (q, uj) = v 2D (q) p\ (q, uj) and 
cxp (— qd) v 2D (q) p\ (q, uj) in sheets 1 and 2, respectively. The resulting potential in 
sheet 2 after screening by the carriers is exp (—qd) v 2D (q) pi (q, uj) / [1 + a (q, uj)], 
which gives rise to an induced carrier distribution in sheet 2, 

„ (q ,„)_ xtaB ,.-^(,) 1 £^L 5r (7) 

In complete analogy, this carrier distribution in sheet 2 gives rise to a carrier dis- 
tribution in sheet 1, 

Pl (q, uj) = X (q, w ) e-^v 2D (q) ft^^tt - (8) 

[1 + a(q,uj)} 

To find the condition for self-sustained fields, normal modes, we let this induced 
carrier density in sheet 1 be the carrier density we started from. This leads to 

2 



1 - e" 2 « d 



«(q,^) 

1 + a(q,w) 



= 0. (9) 



The left hand side of this equation is exactly the argument of the logarithm in 
Eq. ([5]) . That equation was derived using many-body theory. 
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2.2. Results 

The numerical results for the size of the interaction energy between two undoped 
graphene layers in vacuum is shown as the lower straight line in Fig. [5] The bent 
curves are valid for doping concentrations 1 x 10 10 ,1 x 10 n ,l x 10 12 ,andl x 
10 13 cm -2 , respectively, counted from below. The interaction energy is negative, 
leading to an attractive force. The uppermost thick straight line is the classical 
Casimir result for two perfectly reflecting half spaces. We note that the two straight 
lines are parallel and will never cross. The curves representing doped graphene 
sheets on the other hand will, for distances well outside the region covered here, 
approach and join the classical Casimir result. 

3. A Graphene Sheet Above a Substrate 
3.1. Derivations 

This geometry is not well suited for using many-body theory so we go straight to 
an approach using electromagnetic normal modes. 

3.1.1. Normal-mode derivation 

If we have a 2D layer (like a graphene sheet) the distance d above a substrate the 
procedure is very similar to in Sec. (2.1.2). We start with an induced mirror carrier 
density p\ (q, cj), in the substrate. The induced carrier density in the 2D layer is 
given by the expression in Eq. except that now the distance between the mirror 
charge and the 2D layer is 2d instead of d. Eq. ([8]) is then replaced by 




d(A) 



Fig. 2. The attractive interaction energy between two graphene sheets. The lower straight line 
is for undoped sheets, while the four bent curves are for doping densities 1 X 10 10 , 1 X 10 11 , 1 X 
10 12 ,andl X 10 13 cm - 2 , respectively, counted from below. The uppermost thick straight line is 
the classical Casimir result for two perfectly reflecting half spaces. 
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Pi (q,w) = -p 2 (q,w) 



and the condition for normal modes becomes 



£ s H - 1 
e s M + 1 ' 



x _ c -2 q d «(q^) £ s H-i = 



1 + a (q,w) e s (w) + 1 



resulting in the energy 

2 OO ( 

E c (d) = hf^J^ln{l 



^-2qd 



a'(q,uj) e s /(co) — l 
l-\-a' (q.uj) e s f(uj)-\-l 



(10) 



(ii) 



(12) 



3.2. Results 

The result for a graphene sheet above a gold substrate is shown in Fig. [3] The 
dielectric function of gold along the imaginary frequency axis was obtained from 
experimental data extrapolated in a way described in Ref. [T7] and by the use of a 
modified Kramers Kronig dispersion relation (see Eq. (6.75) in Ref. [TCI) 



4. Predicted Power Laws and Actual Distance Dependence 

One way to find fast results for the van der Waals and Casimir interactions betw een 
objects of various shapes is to sum over pair interactions. According to Langbein^SI 
one finds the correct separation dependence but the overall strength is not always 
right. For two parallel thin films this gives the power law dr A (d~ 5 ) in the non- 
retarded (retarded) limit. We found in Ref. [15] that this is not true for a pair of 
2D metallic sheets. We found a fractional power law, d~ 5 / 2 , in the non- retarded 




Fig. 3. The attractive interaction energy between a graphene sheet and a gold substrate. The 
lowest curve is for an undoped sheet, while the the following four curves are for doping densities 
1 X 10 10 ,1 X 10 11 ,! X 10 12 ,andl X 10 13 cm -2 , respectively, counted from below. The uppermost 
thick straight line is the classical Casimir result for two perfectly reflecting half spaces. 
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limit (see the straight part of the solid curves in Fig.[4j) and the power law d~ 3 
in the retarded. It is also interesting to note that the interaction weakens in the 
small separation limit for separations of the order of the Thomas-Fermi screening 
length. This effect is a direct result of spatial dispersion. The importance of spatial 
dispersion in this geometry was also found by GerlacbA^. Later we studied 3D gold 
films of different thickness. We found^Uhat there is a region around the cross-over 




Fig. 4. The attractive interaction energy between two 2D metallic sheets. The curves are for 
doping densities 1 X 10 10 , 1 X 10 11 , 1 X 10 12 , andl X 10 13 cm -2 , respectively, counted from below. 
For large separations the curves follow the classical Casimir result for two perfectly reflecting half 
spaces; for intermediate separations, in the van der Waals region, the curves varies as d~ 5 / 2 ; for 
separations smaller than the Thomas-Fermi screening length the interaction shows a weakening 
due to spatial dispersion effects See Ref. 15 for details. 
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Fig. 5. The attractive interaction energy between two thin gold films. For large separations the 
curves follow the classical Casimir result for two perfectly reflecting half spaces; for small separa- 
tions they follow the van der Waals result for two gold half spaces; in a region in between these 
regimes thin films behave as strictly 2D metal films in the van der Waals regime with a rf — 5 / 2 
separation dependence. See Ref. 20 for details. 
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point between van der Waals interaction and Casimir interaction where thin metal 
films behave as strictly 2D films (See Fig. [5]) Here, in the present work we have a 
system with yet another separation dependence. For a pair of undoped graphene 
layer s th e non- retarded interaction varies as d~ 3 , verifying the prediction by Dobson 

ini I 

et al. 1 ^, and the same power law holds in the retarded regime. This is the same 
power law as for the retarded interaction between two half spaces. The origin of 
the half-integer behavior for the 2D metal sheets is the square root dependence of 
the dispersion curve for 2D plasmons. In doped graphene the plasmon dispersion 
curve attains the square root dependence and for larger separations the interaction 
varies as g? -5 / 2 (see the rightmost part of Fig. [2]). For a graphene sheet above a 
metal substrate we find a more complicated behavior, i.e., no simple power law. 
Summation over pair interactions suggests a d~ 3 dependence in the non-retarded 
regime. 



Table 1. Predictions a la Langbein (within parentheses) and outcome 



Geometry 


van der Waals 


Casimir 


2D metal — metal half space 


d -5/2(a) (d -3) 


rf -3(a) ( d -4) 


2D metal — 2D metal 


d -5/2(a,b) ( d -4) 


d -3(a,b) (d -5) 


graphene — metal half space 


No pure power law (d~ 3 ) 


d~ 3 (d- 4 ) 


graphene — graphene 


d -3 Ccrf) (d- 4 ) 


d~ 3 (d~ 5 ) 


doped 


d -5/20) (cr 4 } 


d-' A (d~ 5 ) 


(a) Bostrom and Sernelius, 

(b) Sernelius and Bjork, 

(c) Bo E. Sernelius, 

(d) Predicted by Dobson et al., 


Phys. Rev. B 61, 2204 (2000). 
Phys. Rev. B 57, 6592 (1998). 
EPL 95, 57003 (2011), (present work). 
Phys. Rev. Lett. 96, 073201 (2006). 





Our calculations were done for T = 0. We have so far refrained from taking finite 
temperatures into account since then one has to use a finite temperature dielectric 
function for graphene, a function that has to be calculated numerically. In Ref.[22] 
one recently used a quantum-field-theory model to derive the retarded interaction 
at finite temperature. One put the main emphasis on the high temperature limit 
(or the large separation limit at finite temperature.) One found that the results for 
two freestanding graphene sheets and for a graphene sheet above a metal substrate 
both agree with the result for two Drude type metal half spaces; this interaction is 
half the strength of the interaction between two ideal metal half spaces and varies 
as dr 2 . The result agrees with that of Gomez-Santos^ but there are some disagree- 
ments as to when this high-temperature result sets in. Another even more recent 
publication^^ treats the retarded interaction between two freestanding graphene 
sheets, undoped and doped, at finite temperature. In the undoped case they find 
for small separations a d~ 3 dependence in agreement wi th o ur results as well as 
with those of Gomez-Santos^ and those of Dobson et aL^D at large separations 
they find a e?~ 4 dependence in disagreement with Gomez-Santos^. In the doped 
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case they find for small separations a d~ 5 ^ 2 dependence and for large separations a 
dr 2 dependence. This is in disagreement with our results in that we find the dr^l 2 
dependence sets in at a much larger separation. We believe the disagreements are 
caused by their treatment of the graphene sheets. They use the formalism for layers 
of finite width. They let each graphene layer keep a small but finite width and let 
the bulk dielectric function be represented by the 2D graphene dielectric function. 
We believe this is not the correct way to proceed. It produces an effective 2D dielec- 
tric function, different from the graphene dielectric function. In Ref.[20]we used the 
bulk dielectric function for gold films of finite thickness. When the film thickness 
goes towards zero the result approaches, as it should, that for a strictly 2D metal 
film. Had we instead used the dielectric function of a 2D metal film for the bulk 
dielectric function the result would have been quite different. 

5. Summary and Conclusions 

We have derived the non-retarded Casimir interaction between two parallel 2D 
sheets and performed numerical calculations for two undoped and doped graphene 
sheets. Spatial dispersion was fully included. Spatial dispersion has strong effect on 
the contribution from the doping carriers but negligible effect on the contribution 
from undoped graphene. We have furthermore derived the non-retarded Casimir 
interaction between a 2D sheet and a substrate. Numerical results are presented for 
an undoped and doped graphene sheet above a gold substrate. In both geometries 
the power law for the separation dependence of the interaction is quite different 
from that predicted by Langbein. The distance dependence is summarized in Table 

m 
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